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We describe a semi- analytical approach to non-linear diffusive shock accelera- 
rrl , tion in the case in which nuclei other than protons are also accelerated. The 

i-H • structure of the shock is determined by the complex interplay of all nuclei, and 

in turn this shock structure determines the spectra of all components. The mag- 
netic field amplification upstream is described as due to streaming instability of 
all nuclear species. The amplified magnetic field is then taken into account for 
O ' its dynamical feedback on the shock structure as well as in terms of the induced 

modification of the velocity of the scattering centers that enters the particle 
^ ■ transport equation. The spectra of accelerated particles are steep enough to be 

compared with observed cosmic ray spectra only if the magnetic field is suffi- 
ciently amplified and the scattering centers have high speed in the frame of the 
background plasma. We discuss the implications of this generalized approach 
1^ ' on the structure of the knee in the all-particle cosmic ray spectrum, which we 

^SJ I interpret as due to an increasingly heavier chemical composition above lO^^eV. 

0^ ' The effects of a non trivial chemical composition at the sources on the gamma 

ray emission from a supernova remnant when gamma rays are of hadronic origin 
[~~. ■ are also discussed. 

^—v ' Key words: supernova remnants; shock waves; Galactic cosmic rays; nuclei; 



knee 



1. Introduction 

A satisfactory understanding of the origin of cosmic rays (CRs) must deal 
with the issue of chemical composition. Probably the most striking instance of 
the role played by the different chemicals is the origin of the knee in the all- 
particle CR spectrum. In a scenario in which the maximum energy of accelerated 
particles scales with the charge of the particles involved, a knee arises naturally 
as a superposition of spectra of chemicals with different charges Ze. Even more 
important, the change of spectral slope on the two sides of the knee is determined 
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by the relative abundance of different chemicals as a function of Z, convolved 
with the effects of rigidity dependent propagation in the Galaxy. 

While there has been much work on the propagation of nuclei in the Galaxy 
(mainly because of the importance it has for the prediction of travel time of CRs 
in the Galaxy and ratios of secondary to primary fluxes), not much attention 
has been devoted to the acceleration of nuclei in the sources. CR acceleration 
in SNRs is believed to take place through the mechanism of diffusive shock 
acceleration, in its non-linear version that allows one to take into account the 
reaction of accelerated particles on the plasma and on local magnetic fields. 
Several versions of the non-linear theory of diffusive acceleration at shocks have 



been developed (see [3J] for a review) , but most of them include only protons 
as accelerated particles. Two noticeable exceptions are represented by the work 
of [lOJ and that of [38| . In both papers the calculations consist of a numerical 
solution of the coupled equations of CR transport and conservation of mass, 
momentum and energy flux of the overall CR plus background plasma. In the 
paper by Berezhko and Volk [lOJ the calculations were not illustrated in detail, 
and it is difficult for us to appreciate the assumptions that were adopted there. 
The spectra from individual sources was found to be very flat (flatter than i?"^ 
at high energies), so that the observed CR spectrum could be recovered only 
by assuming a Galactic diffusion coefficient as steep as D{E) ex E^-"^^, which 
however is not consistent with measurements of the CR anisotropy at the Earth. 



In the paper by Ptuskin et al. 38|, some more details were provided, and the 
authors discussed the important role of the velocity of the scattering centers on 
the shape of the spectrum of accelerated particles, which, as a consequence, is 
here much steeper than that found by Berezhko and Volk [lO]. 

A fully non-linear theory which includes nuclei is made rather complex by 
at least two issues: (1) nuclei change the structure of the shock, making the 
problem harder to tackle and (2) the injection of nuclei in the accelerator is 
more challenging to be modelled than it is for protons, especially because nuclei 
can be produced as a result of dust sputtering. In a non-linear theory, the 
second point clearly feeds back onto the first one. 

It is important to recall that the injection and the non-linear acceleration of 
nuclei in the Earth's Bow Shock (EBS) has been successfully described using a 



Monte Carlo simulation in the pioneering paper in Ref. [25|. These calculations 
showed how the case of the EBS is and remains, 20 years later, the most clear cut 
instance of occurrence of Non-Linear Diffusive Shock Acceleration (NLDSA) in 
collisionless shocks. However some caution should be adopted in extending these 
results to the case of SNRs, in that the physical conditions at the EBS might 
be somewhat different: the lack of electron injection, the fact that magnetic 
field amplification in the EBS appears to remain in the quasi-linear regime and 
the crucial role played by dust sputtering for ion injection in SNRs might be 
considered as possible evidences of such differences. 

The last point is particularly relevant to our problem in that, in order to ac- 
count for the observed discrepancy between the chemical composition of typical 
interstellar medium and CRs in our Galaxy, refractory elements (such as Mg, 
Al, Si and Fe) have to be injected in a preferential way in the acceleration pro- 



cess with respect to volatile elements [35|, |36| . Refractory elements are usually 



trapped in dust grains, and their preferential injection has been interpreted as a 
consequence of the sputtering of the grains when they are swept up by the SNR 
shock. In particular, the processes which lead to the injection of suprathermal 
iron nuclei as a consequence of the sputtering of accelerated dust grains have 
been put forward quantitatively for the first time in j23| . 

Nevertheless, a self-consistent description of this process would require both 
a detailed knowledge of the dust chemistry close to a SNR shock and an accurate 
physical description of the grain sputtering, along with a time-dependent treat- 
ment of the ionization of dust and atoms during their acceleration. The intrinsic 
complexity of such a phenomenon led us (as well as the authors of the previous 
work mentioned above 10|j|38|) to a simplified treatment of the nuclei injection: 



the measured spectra of relativistic particles at Earth are reproduced without 
explicitly taking into account the microphysics either of the nuclei ionization or 
of the dust sputtering (see ^ . 

One might argue that the relatively low abundances of nuclei in the cosmic 
radiation observed at Earth might make their influence on the shock structure 
negligible, so that in describing the acceleration process, the shock structure 
could be treated as determined by protons alone, with nuclei behaving as test 
particles. However, after correcting for propagation effects, it is easy to show 
that the nuclear contribution to the total pressure and magnetic field amplifica- 
tion in the vicinity of a typical supernova remnant shock may be as important 
as that of protons. 

In this paper we describe the generalization of the non-linear theory of DSA 
developed by Amato and Blasi [71, |8| and Caprioli et al. j21| to include nuclei 
of different charges. We calculate the spectrum of all species as accelerated at 
the shock and the structure of the shock (including magnetic effects) induced 
by them. We also comment on the implications of acceleration of nuclei on 
the spectra of secondary products of particle interactions, especially gamma 
rays. Finally we show the all-particle CR spectrum at Earth resulting from this 
calculation. 

The calculations discussed here are semi-analytical, and from the computa- 
tional point of view very inexpensive. This allows us to explore a wide region 
in parameter space, which is particularly important when dealing with the goal 
of explaining the CR spectra and chemical abundances observed at Earth. For 
simplicity we consider here only supernovae exploding in a homogeneous inter- 
stellar medium (ISM). While more complex situations can be treated in the 
context of our formalism, they introduce a wide range of new and hardly acces- 
sible parameters, which overshadow the main physical results. We will comment 
further on this point whenever we deem it necessary. 

The paper is organized as follows: in §[5]we illustrate the generalization of the 
equations and solution techniques following the non-linear theory of 0, H, l2l| . 
In §[3] we comment on the abundance of nuclei in CRs and at the sources, and 
how it relates to the injection of the nuclear component at the shock. In §|4]we 
illustrate our results for the spectra of accelerated particles and the structure of 
the shock. In §[5] we discuss the implications of the presence of accelerated nuclei 



at the shock for the prediction of a gamma ray flux, as due to production and 
decay of neutral pions. In §[6] we compare our findings with the CR flux detected 
at Earth, discussing the spectra of individual chemicals and the appearance of 
the knee a.t E ^ 10^ GeV. We discuss our general results and we compare them 
with previous findings in § [T] 

2. Equations and Solution techniques 

In this section we generalize the semi-analytical formalism developed in |7|, |8|, 



2lj to the case in which nuclei heavier than Hydrogen (hereafter simply Heavy 
Nuclei, HN) are also injected and accelerated at a stationary, plane, parallel 
(background magnetic field parallel to the shock normal) , newtonian shock wave. 
We label with a subscript i quantities referring to different chemical elements, so 
that the convection-diffusion equation for the isotropic part of the distribution 
function, fi{x,p), reads, for each species [see e.g. 39]: 

,dft{x,p) 



dx dx 



D,{x,p)- 



dx 



l^^+*<-^ m 



Here Di{x,p) is the parallel diffusion coefficient, which may depend both on 
space and momentum, u{x) = u{x) + vw is the total velocity of the scattering 
centers in the shock frame, given by the sum of the fiuid velocity u{x) and 
wave velocity wvy, and Qi{x,p) is the injection term. The shock is at a; = 
and subscripts 0, 1 and 2 label quantities taken, respectively, at the upstream 
free-escape boundary x — xq, immediately upstream and downstream of the 
subshock. 

An especially important issue, when taking into account nuclei, is that of 
particle injection, as stated in the introduction and as should become clear be- 
low. Unfortunately, the microphysics of this process is not yet fully understood 
even for the case of protons alone, and much more so for nuclei. In the follow- 
ing we simply assume that protons are injected from downstream via thermal 
leakage as described in [ij], while the injection of HN is tuned in such a way as 
to reproduce the relative abundances observed in the CRs. We do not account 
for the details of the HN injection, likely related to the complex physics of the 
dust sputtering process. Indeed, the required preferential injection of HN is 
likely related to the fact that partially ionized heavy particles (i.e., thermalized 
particles with large mass/charge ratios) have large Larmor radii and are hence 



preferentially injected in a thermal leakage scenario (see [2J|), and/or in the fact 
that refractory nuclei can be efficiently injected via dust grain sputtering [23| 
(we will discuss these points in S|3|). 

More precisely, we assume that all protons with momentum p > Pinj,ii have a 
large enough Larmor radius to cross the shock (subshock) from downstream and 
start being accelerated. Since the shock thickness is expected to be of the order 
of the Larmor radius of particles with thermal momentum pth = v^^tihAbTh^ 
(where rH,2 is the downstream proton temperature and ks is the Boltzmann 
constant), we take Pinj,}i — Ch Pt/i,H, with ^ '-^ 3 — 4. Furthermore, at a given 



momentum p the Larmor radius of a HN with charge Z^e is a factor Zi smaller 
than that of a proton. Hence, in this scheme, it is very natural to assume that 

Pinj^i — ^i Pinj.H- 

Clearly, if nuclei are not completely ionized or if they are injected via dust 
sputtering this simple recipe may fail to describe the low-energy tail of the CR 
distribution. Nevertheless, since the spectra of the relativistic particles observed 
at Earth are consistent with power-laws extending from a few GeV/nucleon up to 
the knee, it seems reasonable to assume that injection always occurs at energies 
between thermal and mildly relativistic. Therefore, the spectra we work out 
are expected to be accurate throughout the entire energy region relevant for the 
SNR emission, and for the measured Galactic CRs as well. 

The injection term in Eq. [T]can be written as 



.,{x,p) = Q,^i{p)6{x) 



Airp: 



Sip-pinj.i)S{x), 



(2) 



inj,i 



where rji has the usual meaning, i.e. the fraction of particles of each species 
crossing the shock that is injected in the acceleration process. We will discuss 
in 331 how to determine reasonable values of 77^ as inferred from what is measured 
at the Earth. 

We solve Eq.[l]along with the spatial boundary condition fi{xo,p) — 0, which 
mimics the presence of a free escape boundary upstream placed at x = a;o- For 
each chemical element, the distribution function fi{x,p) and the escape flux 
(j>esc,i{p) can be written as [seel21|: 



f^{x,p) = /s/,,j(p)exp 



, u{x') 
D^{x',p) 



1 



W,{x,p) 



W,fl{p) 
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Dt{x,p) — 



uofsh,i{p) 



where we have introduced the functions fsh,i{p) — /i(0,p), and 



Wi{x,p) =Mo 
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(4) 



(5) 



The distribution function at the shock reads: 



fshAp) 



Vino 



3Rf 



inp^ 



RtotUp,i{p) - 1 



exp ■ 
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(6) 



where we used 



UpAp) = Ui 



dx 



dU{x) fr{x,p) 



(7) 



'xo do; fshAp) 

and introduced the subshock Rsub = ui/u2 and the total Rtot = uq/u2 compres- 
sion ratios for the fluid. In this approach in which the waves move with velocity 
VA with respect to the fluid, the cosmic rays which scatter against them feel an 



effective velocity given hy u + va and hence a subshock and a total compression 
ratios given by 



R. 



U2 + VA,2 '' 



Rt, 



Up + VA,0 
"2 + VA,2 ' 



(8) 



As we discuss below, the particles may feel compression ratios larger or smaller 
than the fluid, depending on the relative sign of u and va- 

The convection-diS'usion equation is coupled with the standard conservation 
equations for mass, momentum and energy flux, where the cosmic ray terms are 
intended to be summed over all the elements. In particular, the conservation of 
momentum, after dividing by the bulk pressure PqUq, reads: 



U{x) + Pg{x) + Pc{x) + P^{x)^l + 



iM^ 



(9) 



where the gas pressure, if the heating in the precursor is purely adiabatic, can 
be written as 

If turbulent heating is also at work converting a fraction C^ of the magnetic 
pressure into thermal energy, the pressure of the gas would instead read [see 
§4.2 of[l3: 



p,{^)^^^^[^ + ':h{x)] 



iJ(x)- 7(7-1) 



Ml_ 
Ma.o 



1- [7-^+1/2 (x) 



7 + 1/2 



The pressure in cosmic rays is, as usual. 



A P+OO 

Pe(x) = J2 P^(^) = ^ E / '^PP'' "(P) ^^(^'P)- 



(11) 



(12) 



Finally the pressure in the form of magnetic turbulence, which has been shown 
to play a key role in the shock dynamics [16[, can be approximated, assuming 
that only standard Alfven waves are generated via resonant streaming instability 
and for Mq, Ma > 1, as [III 



pux)^{i-ou{xr^/' 



1 - U{xf 
4Ma,o 



(13) 



A solution of the problem of NLDSA can be obtained through the same recursive 
method described in [21|. We summarize it here. 

• Start from a guess for U{x) and for the distribution function fi{x) {e.g. the 
test-particle solution) and fix a value of U* = Rsub/ Rtot- This corresponds 
to a value of P*i, through Eq. |9l since Pg and P^ depend only on U. 

• Compute the new distribution functions fi{x,p) by using U{x) and fi{x,p) 
from the previous step in Eqsl5][7] 



• Compute the CR pressure through Eq. [T^l In general, the value of P^i 
that is obtained will be different from P*i. Then Pc{x) and fi{x,p) are 
renormalized by multiplying them by a factor A = P*-^^/Pc.i- 

• Compute and update the fluid velocity profile U{x) using the new Pc{x) 
in Eq. M 

• Iterate the steps above until A does not change between a step and the 
following. 

Starting from an arbitrary [/*, in general one has A 7^ 1 after convergence is 
reached. The process must then be restarted with a new choice of U* until 
convergence to A = 1 is achieved (within a prescribed level of accuracy). The 
fluid profile and the distribution functions so obtained are, at the same time, 
solutions of both the convection-diffusion and the conservation equations, i.e. 
solutions of the full problem. 

It is worth recalling that from the computational point of view this approach 
to non-linear shock acceleration is the most efficient presented so far in the 
literature: for a given set of parameters it takes less than one minute (on a 
laptop) to calculate the distribution function fi{x,p) of all accelerated particles 
and the space dependent structure of all thermodynamical quantities in the 
shock region. 

3. Abundances of heavy elements 

The most challenging aspect of a non-linear calculation of diffusive shock 
acceleration in the presence of accelerated nuclei is the tuning of the relative 
abundances at the source in a way that may fit the abundances observed at 
the Earth after propagation. The difficulty is manyfold: first, the non-linearity 
of the problem makes it very difficult to establish once and for all what the 
required relative abundances are; second, the spectrum observed at Earth from 
an individual SNR is the superposition of the instantaneous spectra of particles 
escaping the remnant at different times, and the relative abundances are, most 
likely, time dependent. One final complication to keep in mind in establishing a 
connection between source and local abundances of ions comes from the spalla- 
tion processes suffered by nuclei during propagation in the Galaxy. This simply 
means that the relative abundances of chemicals in a given remnant at a given 
time cannot be univocally inferred from the CR abundances that we measure 
at Earth. 

It is worth recalling that the situation is rather different for the EBS, where 
real-time measurements of the spectra of different elements are available, thus 
allowing a more detailed interpretation of the ongoing physical processes. In 
this case, in fact, the injection of both protons and nuclei via thermal leak- 
age provided a successful explanation of the EBS properties as due to efficient 
NLDSA ^2^. 

On the other hand, the same thermal leakage scheme seems unsatisfactory 
when applied to the case of SNRs as the sources of Galactic CRs. In fact, a 



thorough study of the elemental and isotopic ratios in CRs compared with solar 
ones shows that refractory elements are preferentially accelerated with respect 
to volatile ones. We refer to the original works 35|, |36j (and to a recent review 



41|) for a comprehensive discussion on these topics. 

Here we only want to stress that this fact has been explained by invoking 
the fundamental role of the interstellar dust in the injection of nuclei, as put 



forward in 23|. In few words, when a grain of refractory material is caught by a 
SNR shock, it is efficiently accelerated by crossing the shock many times thanks 
to its large Larmor radius. When its velocity reaches about 0.01c it begins to 
sputter and to release atoms. These have the same velocity of the parent grain: 
they are hence well suprathermal and therefore can be easily injected and take 
part in the acceleration process. Such an effect can also account for the fact that 
among refractory nuclei there is no preferential acceleration of heavy elements 
with respect to lighter ones, which is instead the case for volatile elements. 
This latter evidence is usually explained as a consequence of the fact that heavy 
elements may be ionized only partially and hence preferentially injected because 



of their relatively large rigidity (see [22|, |2J]) 



In this section we illustrate for simplicity how to estimate the relative nor- 
malization between source spectra of different ions in a test-particle approach. 
Such an estimate is expected to hold for energies above a few GeV/nucleon, 
where the spectra detected at Earth are very close to power-laws. The situa- 
tion may be more complicated in the region between slightly suprathermal and 
weakly relativistic particles because of the effects of partial ionization and dust 
sputtering discussed above. 

The resulting abundances, not too dissimilar from the exact ones derived in 
^ are used to make several important points, mainly concerning the instan- 
taneous spectra of relativistic particles (for instance to address the dynamical 
role of nuclei at the shock, and their effect on the generation of secondary radi- 
ation). This discussion can be carried out regardless of a detailed knowledge of 
the non-relativistic region of the spectra. 

The starting point of our estimate is represented by the observed relative 
abundances in CRs at Earth [see e.g.|30|, for a recent review of the observational 
results] . In a test-particle approach the spectrum of particles accelerated at the 
shock can be written as a power law with slope /3: 

Crji f P 



'•"'•sattei ■ <"' 



with C the same constant for all i. Let us also assume Pinj,i — ca Pinj,ii- In the 
following we shall assume that the injection is charge dependent, with ai = Zi, 
being Z^e the electric charge of ions of type i. The number of particles with 
momentum > p can be estimated as Fj(> p) w Anp'^fi{p) — C?]i{p/pinj^i)^~^ ■ 
The spectrum at Earth is affected by propagation as 



\ 3-/3 

n^ (p) « F, {p)n (p) = Cm ( -^ Zfm (p) , (15) 

^Ptnj,i / 



Element Flux(l TeV) 7^ K~^ ^ 

H (8.73 ±0.07) X 10-^ 2.71 ± 0.02 1 1 

He (5.71 ± 0.09) X 10-2 2.64 ±0.02 0.90 0.28 

C-N-0(*) - 5.4 X 10-2 (*) ^ 2.64(*) 0.85 0.03 

Mg-Al-Si ~ 1.7 X 10-2 ^ 2.66 0.25 3.2 x 10-^ 

Fe (2.04 ± 0.26) X 10-2 2.59 ± 0.06 0.41 1.6x10' 



3 



Table 1 : Flux at 1 TeV (in units of [ni^ sr s TeV] - ^ ) , spectral slope (7^ ) and ratio relative to 
H at 10^ GeV (Km) measured at Earth for the most abundant elements in galactic CRs. CR 
data are from [2911 . table 7, except the ones denoted with (*) which are from [3^. C, N and 
O are taken as one effective element with Zcno = 7 and j4cno = 14- The same is done for 
Mg, Al and Si, with ^MgAlSi = 13 and AjvigAlSi = 27. 

where we used the fact that the confinement time r^ of nuclei of charge Zi is Zf 
times larger than that of protons at given momentum. This reflects the rigidity 
dependent nature of particle propagation in the Galaxy, namely the fact that 
the diffusion coefficient is D{p) (x R^ — {p/ZiY [R is the particle rigidity). The 
simple scaling in Eq. [12] holds as long as spallation in the ISM can be neglected, 
therefore it should only be used at sufficiently high energies (we use p* = 10^ 
GeV/c as normalization point) so that the escape time from the Galaxy is much 
shorter than the spallation time scale for all the relevant species. 

The ratio of abundances between ions and protons at the same momentum 
p* as measured at the Earth is therefore: 

i^^H = ^ = ^Z,f f^^^^V"^ = ^Zf+^-^ (16) 

In this expression, the ratio KiY^ comes from the measured spectra and, at 
least in the context of this simple (leaky-box-like) approach, /3 + 5 is fixed by 
observations to be ^ 4.7. The ratio of the Ty's to be used in the calculations can 
therefore be easily estimated from measurements. The values resulting from 
this procedure are illustrated in Tab. [T] In the same table we also show the 
flux of different chemicals at 1 TeV and the spectral slope at Earth as inferred 
using the so-called Poly-Gonato parametrization by Horandel [29| . These slopes 
reflect the presence of spallation and most likely the superposition of different 
types of sources (see discussion in iJTl) . 

As stressed above, the procedure just outlined is based upon a test-particle 
picture and should only be used as a qualitative estimate of the efficiency of 
acceleration of the different chemicals. Many factors can affect the result of this 
simple estimate: (1) the non-linear diffusive acceleration at shocks does not lead 
to power law spectra; (2) the values of rji can be (and in general are expected 
to be) time dependent, which is especially relevant since the spectrum of CRs 
observed at Earth results from the superposition of instantaneous spectra of 
particles escaping the SNR at different times; (3) there are different types of 
SNRs exploding in different environments, which leads to different ratios at the 
sources; (4) spallation changes the spectra in a complex way, and the difference 
in the spectral slopes of different chemicals as reported in the table above can 



partially be due to these reactions, although it is likely that the differences also 
reflect the contributions from SNRs in different ambient media; (5) the diffusion 
coefficient is not strongly constrained by the secondary to primary ratios and 
the slope S could well be in the range 0.3 — 0.7, though the upper end of this 



range leads to problems with anisotropy (e.g. [28[); (6) although the spectrum 
observed at Earth is very close to a power law, different SNRs or the same 
remnant at different stages of its evolution may have both different spectra and 
different values of the relative abundances. 

In 21 and fj5] we adopt the simple procedure outlined above to estimate the 
abundances of ions in a source at a given time, while in SJ51 where we calculate the 
spectrum of the different chemicals at the Earth and the all-particle spectrum, 
we follow a somewhat more refined procedure. It is however wise to keep in 
mind that, for all the reasons mentioned above, it is extremely difficult, if at all 
possible, to infer in a realistic way the abundances at the sources that may fit 
the whole set of data available at Earth. 



4. Particle spectra and shock hydrodynamics 

We consider a shock moving with velocity uq = 4000 km/s in a homogeneous 
medium with particle density 0.01 cni"-^, temperature Tq = lO^K (sonic Mach 
number '^ 34) and magnetic field Bq ~ 5/iG aligned with the shock normal. This 
choice of parameters corresponds to a 2000 year old SNR with radius Rsh ^ 14.4 
pc, i.e. a SNR at the beginning of its Sedov- Taylor stage for a SN explosion of 
10^"'^ erg and an ejecta mass of 1.4 solar masses. It is worth recalling that the 
highest cosmic ray energy is thought to be achieved at this evolutionary stage 

The free-escape boundary is placed at xq ~ Q.2Rsh upstream of the shock 
and the diffusion coefficient is taken as Bohm-like, 

D,{x,p)^\v{p)-^ (17) 

3 ZiB{x) 

in the amplified magnetic field at the shock position, namely B{x) = Bi = 
•\/87r/9oWgPiu,i upstream and B{x) = B2 — RsubBi downstream. 

In Fig. [1] we show the spectra of accelerated particles and the the quan- 
tities related to shock hydrodynamics, obtained through the iterative method 
described in Sj2l in a case of efficient particle acceleration (we used ^h = 3.8, cor- 
responding to ?7h = 5.7 X 10~^ in Eq. [2|). Notice that the gas pressure upstream 
is very low and lies outside the plot boundaries. 

The most noticeable feature is the fact that, for the standard abundances 
deduced in ^ the dynamical role of nuclei overall is twice as important as 
that of protons: at the shock position the pressure of accelerated protons is 
Pjj — 0.05, in units of the ram pressure far upstream, while the pressure in 
the form of relativistic HN is about 0.1 (right panel of Fig. [T]). The latter is 
mainly associated to He nuclei (Pro — 0.045) but also CNO and Fe nuclei carry 
non-negligible fractions of the momentum in the form of accelerated particles 
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Figure 1: Left panel: all-particle spectrum (thick line) and spectra of individual elements. 
Right panel: spatial dependence of the hydrodynamical quantities (the gas pressure upstream 
is very low and lies outside the plot boundaries). 



(-PcNO — 0.017 and Pfc — 0.003). The relative importance of HN with respect 
to protons should not be affected by the exact shape of the spectra between 
Pinj,i/Zi and a few GeV in rigidity (which is related to the injection details 
discussed above), provided that the spectra are not much steeper than p~'^ 
and/or that the maximum momentum is larger than few GeV/c. These peculiar 
conditions might only be met rather late in the remnant evolution, and definitely 
not as early as at the beginning of the Sedov stage, i.e., when the acceleration 
efficiency is at its maximum. 

For comparison, in the right panel of Fig. [T] we also show the pressure asso- 
ciated to the thermal plasma and to the magnetic turbulence (taken as homo- 
geneous upstream) and the normalized fluid velocity as well. The latter shows 
the spatial profile typical of cosmic ray modified shocks, with a precursor, a 
weak subshock with compression ratio Rsub — 3.9 and a total compression ratio 
larger than 4, namely Rtot — 4.5. 

The resulting spectra of the accelerated particles reflect the modified dy- 
namics, being steeper than the test-particle prediction, ~ p~^, at low energies, 
and flatter than p~'^ at the highest energies. The finite shock size (spatial 
boundary), which allows for particle escape from upstream, induces a rigidity- 
dependent cut-off in the spectra of the various species, leading to maximum 
momenta scaling with atomic charge (left panel of Fig. [1} . More precisely, since 
the diffusion is rigidity-dependent, and since in diffusive shock acceleration the 

spectral slope at a given momentum, qi{p) = °sjt(P) ^ depends only on the 

compression ratio actually experienced by the particles with that momentum, 
we have the simple scaling: qi{p/Zi) — qaip)- In principle, when the relative 
normalization and the ionization state are fixed, and when the CR pressure is 
dominated by relativistic particles, one could use this simple scaling to calcu- 
late the distribution functions of all the elements at the shock, their spatial 
dependence via Eq. |3]and finally the total pressure in cosmic rays Pc{x) and the 
whole shock structure. However, within the fast semi-analytic formalism out- 
lined in 321 the computational effort is far from being a severe issue demanding 
numerical optimization. 
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For comparison, we solved the set of equations describing a shock with ex- 
actly the same characteristics, except for the fact that no HN were considered. 
We found Pn ~ 0.07 at the shock. Comparing this result with those in the right 
panel of Fig. [U the conclusion is that when acceleration of HN is taken into 
account, the pressure in accelerated protons is decreased, but the total pressure 
in accelerated particles is larger, and the shock correspondingly more modified. 
Indeed, in the absence of HN, the compression ratios turn out to be Rsub — 3.95 
and Rtot — 4.2. At the same time, also the magnetic field amplification via 
resonant streaming instability of the accelerated particles gives different results 
when the acceleration of HN is included. We find that the proton-only case re- 
turns a downstream magnetic field of ~ 33/iG, while with HN B2 — 47/zG. This 
result is easily understood given the dependence of the amplified magnetic field 
on the pressure of accelerated particles. All these effects contribute to make the 
case for the necessity of taking HN into account properly: the acceleration of 
HN cannot be treated in a test-particle approximation, nor can it be linearly 
added to a proton only case. In order to correctly describe the shock, all species 
need to be taken into account in a fully non-linear calculation. 

A crucial point to stress is that the mild shock modification and the cor- 
respondingly steeper spectra found here and illustrated in Fig. [T]are the con- 
sequence of the assumption of scattering centers moving with a wave velocity 
that equals the Alfven velocity calculated in the amplified magnetic field. As 
discussed by Caprioli et al. [20| , without this assumption the spectra produced 
in SNRs are very fiat and can hardly be related to the CR spectrum observed 
at Earth, unless, as done by Berezhko and Volk [lOJ, one assumes a very strong 
dependence of the diffusion coefficient in the Galaxy on energy, D{E) oc E'^-'^^, 
which however leads to too large anisotropy at the knee compared with obser- 
vations. It is therefore of the highest importance to realize how the main result 
in this type of calculations is actually due to one of the aspects that are least 
known. 

The fact that taking into account the velocity of the scattering centers could 
lead to steeper spectra was already recognized in the early literature on the 
topic [see e.g.[9|, but the issue has become much more discussed in recent times, 
after discovery of highly amplified magnetic fields in the shock region. There 
are two levels of the problem: (1) we are unable to determine the velocity 
of the waves responsible for the scattering from first principles, especially in 
the case of relevance for us in which magnetic field is strongly amplified by 
CR induced instabilities. If one assumes (unrealistically) that even in these 
extreme circumstances the waves remain Alfven waves, namely with magnetic 
field perpendicular to the (much smaller) background magnetic field, then the 
Alfven speed is well defined, but its value is so small to induce negligible effects 
on the spectrum of accelerated particles. On the other hand, one might expect 
that an eflFective value of the wave speed is close to the Alfven speed calculated 
using the strength of the amplified field as reference value. In this case the Alfven 
speed is much larger and its effect on the transport equation is not negligible. 
(2) Even if we knew how to calculate the wave velocity (for instance by assuming 
that it corresponds to the Alfven speed in the amplified field), the helicity of 
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the waves is unknown, being related to the type of wave, the mechanism that 
generated them and the transport of the waves through the shock (reflection 
and transmission). In other words, it is not easy to calculate the wave velocity 
with respect to the shock surface. 

Since the slope of the CR spectrum depends on the compression ratio of 
the scattering centers, as felt by the particles, it is clear from Eq. [5]that, when 
vw is not negligible compared with u, the actual spectrum may be flatter or 
steeper than the standard prediction depending upon the relative sign of u and 
Vw, both in the upstream and in the downstream regions. If vw is calculated in 
the amplifled magnetic fleld, it typically turns out to be a not-neglibile fraction 
of u, and thus the hclicity of the waves strongly affects the shape of the CR 
spectrum. 

While upstream the turbulence is very likely generated by the CR gradient, 
and thus vw is expected to have sign opposite to u, downstream the waves may 
retain the same helicity (as one would expect based on the transmission and re- 
flection coefRcients at the subshock discontinuity appropriate for Alfven waves), 
or rather may be fully isotropized {vw = 0). Another possibility, recently ex- 



plored by Ptuskin et al. [38[ , is that downstream waves may also be generated 
via streaming instability by the advected CRs, whose gradient is induced by 
the postshock evolution and in particular by adiabatic losses in the expanding 
shell: this scenario leads these authors to assume that in the downstream vw 
has the same sign as u. This mechanism may be at work in SNRs, but one would 
expect that it should act in the downstream plasma on hydrodynamical spatial 
scales. It is therefore hard to assess the relevance of this process as compared 
with that of the turbulence produced upstream and then advected through the 
subshock in the downstream region. In this work we conservatively assume that 
downstream vw,2 = 0, which leads to a less dramatic dependence of the result 
on the Alfvenic Mach number with respect to the assumptions of [38j, even if 
in both cases Ma is calculated in the amplified magnetic field. 

Moreover, in our model the strength of the amplified magnetic field is not 
fixed a priori but is rather a non-linear output of the full calculation, depending 
on the efficiency of the CR acceleration through Eq. [T21 In this sense our 
approach includes a self-regulating effect which returns almost the same result 
in terms of particle spectra and shock modification for a wide range of values 
of the Alfvenic Mach number at upstream infinity, Ma,o, and of the particle 
injection efficiency £,i. 

5. 7— Rays from accelerated nuclei 

In the typical environment of a SNR, in addition to playing an important 
role on the shock dynamics, accelerated nuclei may also give a contribution in 
terms of 7— ray emission via tt'^ production in nuclear interactions with thermal 
protons. In order to calculate this emission, we model each accelerated nucleus 
with atomic number Ai and energy E as an ensamble of Ai protons, each with 
energy E/Ai , and adopt the parametrization of the nuclear proton-proton cross 
section worked out by Kamae et al. |32l |. 
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Figure 2: Spectra of 7-rays due to the decay of tt s produced in nuclear interactions (thick 
Hne). The partial contributions by each element are also shown. 



A simple estimate of the HN contribution to pion production can be obtained 
as follows. Let us assume that a nucleus with energy Ei produces monochro- 
matic photons with energy E^ such that Ei = x^iE-f, with x ~ 10. The photon 
spectrum produced by nuclei of species i would thus be: 



dE., 



N^'iE,) ex A,N,iE,)—^ = xA^,N,[E,) 



dE^ 



(18) 



For spectra Ni{Ei) oc E'j ' up to an energy Emax,i we have: 






= K,,A^ 



2-q 



(19) 



which means that, for the standard case 9 = 2, all the species contribute to the 
7-ray flux proportionally to their abundances at any photon energy below the 
minimum of Emax,i/ {AiX). Since typically Emax,i oc Zi, the photon spectrum 
of any HN is cut-off at an energy which is a factor 2 lower than for protons. It 
is also interesting to notice that, if the accelerated particle spectra are flatter 
than E~'^, the contribution of HN with respect to protons is boosted by a factor 
A^~'. This turns out to be as large as a factor ^ 5 in the case of Fe nuclei 
accelerated with a spectrum ~ E~^-^. There may be stages in the evolution of 
a SNR when the shock is strongly modified and such hard spectra may appear. 
The spectrum of 7— rays, calculated using the particle spectra shown in 
Fig- El is illustrated in Fig.[5J The HN contribution to the 7— ray flux turns out 
to be dominant over that of protons, with a total predicted flux which is about 
a factor 2.5 larger than the standard prediction for the case when accelerated 
particles are only protons. Moreover, since the contribution of all species of 
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HN is truncated a factor 2 below that obtained with protons alone, the shape 
of the high-energy end of the photon spectrum is somewhat different when HN 
are taken into account. Although this may in principle represent a spectral 
feature flagging for the presence of accelerated nuclei, we find it unlikely to be 
detectable, given the many intrinsic systematics of experimental and theoretical 
nature and given the topology of the emitting regions which in general are 
complex, so that different regions of the remnant (where the highest energy 
could be slightly different) contribute to the flux along the same line of sight. 

For the modeling of the gamma ray emission from SNRs a crucial parameter 
is the gas density in the shock region. Moreover, in the context of non-linear 
theory, the temperature of the downstream plasma is an output of the problem 
and the thermal emission from downstream (continuum due to bremsstrahlung 
and lines from non-equilibrium ionization) can be calculated. The latter, scaling 
with the square of the plasma density, is more sensitive than the gamma ray 
emissivity to the value of the gas density. In some cases the gas density is inferred 
from gamma ray observations, with the strong assumption that the gamma ray 
emission is due to production and decay of neutral pions and checked versus 



the signal in the form of lines of thermal origin [see e.g. |26|, for the case of 
RX J1713. 7-3946]. The results shown in this section illustrate the fact that the 
inclusion of ions in the calculations of the gamma ray emission from a SNR 
leads to an estimate of the ambient density which may easily be a factor ^2 — 4 
smaller than when the gamma ray emission is calculated using only protons. 
The thermal emission decreases correspondingly by ^ 4 — 16. Claims on the 
detectability or non-detectability of thermal emission should take into account 
the intrinsic uncertainty associated with this phenomenon. 

6. The Galactic spectrum of CRs 

In this section we calculate the total spectrum of CRs accelerated by a class 
of benchmark SNRs during the different stages of evolution, as described by the 
analytical evolutionary scheme of [40| . Our reference SNR has the same environ- 
mental parameters as in SJH and we follow its evolution from 0.1 to about 22Tst, 
where Tst — 2000yr is the beginning of the Sedov- Taylor stage. We assume that 
the ambient ISM is homogeneous, hot and rarefied (Tq = 10^ K, no = 0.01) in 
order to mimic the explosion of a core-collapse SN into the hot bubble excavated 
by its pre-SN wind, since about 85% of Galactic SNe are expected to be of this 



type [type II -|- type Ib/c, see e.g. l27[. Such a choice should also be represen- 
tative of the typical hot medium inside superbubbles, where most [about 75%, 
according to |27| of the core collapse SNe are found. We do not account here 
either for the presence of a dense, cold wind, occasionally produced during the 
Red Supergiant stage of the progenitor, or for the finite size of the hot bubble. 
These ingredients add to the problem of particle acceleration numerous difficul- 
ties and require the introduction of several almost unconstrained parameters, 
like for instance the mass-loss rate and the total mass blown in the wind, the 
temperature and the density profile and, more important, the strength and the 
topology of the magnetic field. By the same token, here we did not include the 
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contribution of type la-like SNe, which are expected to be less frequent and to 
explode into the denser and colder homogeneous ISM. 

It is important to stress that our decision to not include these situations 
in our calculations is not due to technical problems of our computation tech- 
nique, but only to our belief that the additional uncertainties that would follow 
overcome the benefits of adding a class of sources to our calculations, thereby 
making the results not more, but in fact less trustworthy. Nevertheless, some of 
these pieces of information may be available when dealing with a given SNR: in 
this respect, our approach could be easily modified in order to include a more 
complex evolution of the remnant accounting also, e.g., for the presence of ei- 
ther a dense Supergiant wind or the final encounter between the forward shock 



and the ordinary ISM (see for instance [37| for an analytical treatment of the 
evolution of different SNR types). 

The spectra of the accelerated particles are calculated as described in [20|, i.e. 
taking into account the instantaneous escape flux from the upstream boundary 



see 



18[ , the advection in the downstream region, which leads to adiabatic losses 
as a consequence of the shell expansion, the escape of particles from downstream 
because of shell breaking and/or inhomogeneities in the circumstellar medium 
[th e fraction of downstream escaping particles is taken as 10%: see also §3.3 of 
I20I for a discussion of this point] . The relative abundances at the sources are 
iteratively adjusted in order to fit the fluxes measured at Earth. 

The damping of the magnetic field is heuristically taken into account by as- 
suming that a fraction ( of the generated turbulence is damped into gas heating 
(Alfven heating), with C(i) = 1 — exp [— Mo(i)/'Udamp], where the velocity Mdamp 
at which the effect becomes relevant may range between about 200 and 1000 
km/s, depending on the details of the turbulence generation and damping [see 
I37I for a study of this topic] . We checked a posteriori that our findings depend 
only weakly on the choice of Udamp- 

The convolution over time of the spectra of accelerated particles injected into 
the Galaxy returns the source spectrum, which has been corrected by accounting 
for the propagation in the Milky Way. We adopt a simple leaky box model of 
the Galaxy, which is taken as a cylinder with radius and half-heigth equal to 10 
and 3.5 kpc respectively, where 3 SN explode each century. The escape time is a 
function of the nucleus rigidity R as inferred by standard secondary-to-primary 
measurements and refiects into a grammage Xescip) — 7.3{p/Z/10GY)~^l3{p) g 
cm^^, with 6 = 0.55 and /3(p) the dimensionless speed of a nucleus of momen- 
tum p. Particle losses due to spallation against interstellar nuclei during the 
propagation is also taken into account and described as in §4 of [3l| . 

The results for the different chemicals are shown in Fig. |3l as plotted on 
top of the experimental data collected in [15|. The agreement is very good, 
and both the spectra and the absolute and relative normalization are consistent 
with the data from direct and indirect experiments able to resolve the chemical 
composition of the CR fiux at Earth. Notice that some discrepancies at low 
energies are due to the fact that the results of our calculations have not been 
corrected for solar modulation. 

It is interesting to notice that the spectra of nuclei at Earth have roughly the 
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Figure 3: CR flux measured at Earth for the different chemicals considered in the text [the 
experimental data are from fig. 9 of|l5l. where the proper references to the sigle experiments 
can be found as well]. The thin solid lines in the H, HE and Fc panels correspond to the 
Poly-Gonato fits |29| . ,see also[l]. 
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Figure 4: All-particle CR flux at Earth in the knee region. The thick line represents our model 
output, while the thin lines guide the eye by showing two power laws, oc E^ ' and oc Eq ' 
below and above the knee, respectively [data from [isl. fig. 8]. The partial contributions from 
H, He, CNO, MgAlSi and Fe are also shown. 



same slope at given rigidity, except for the effect of spallation, which produces 
the low energy turnover. For iron nuclei the flattening due to spallation is visible 
up to energies of ~ 10^ GeV. This conclusion should however be taken with 
caution: the recent data from the balloon-borne experiment CREAM-II [33j 
suggest that spectra of nuclei are somewhat flatter than the proton spectrum. 
This might be a mild signature of the contribution of different classes of sources 
operating in environments with different chemical composition and contributing 
to the CR spectrum in a different way than described here. 

In Fig. |4]the all-particle spectrum in the knee region is shown (thick line), 
again compared with data from [15| . For comparison two power-laws are also 
drawn (thin lines), corresponding to the flducial slopes of the flux below {E^'^ ) 
and above {Eq^'^) the knee. 

The all-particle spectrum that we calculated is, as expected, somewhat 
steeper than the standard best fit to the data above the knee (oc Eq^'^). For 
this effect, also found by other authors, several possible explanations have been 
proposed: in §2 of the review by Hillas [28[ the author discussed the need for 
an "extended tail" in the Galactic flux in order to explain the CR data in the 
region between 10^ and 10^ GeV (see fig. 2 of that paper), and suggested that 
such a high-energy feature might be granted by a proper account of a peculiar 
class of Type II SNe. A similar scenario has been put forward by Ptuskin et al. 
[38| , who successfully explained this region of the spectrum with rare (only 2% 
of all Galactic SNe) but very energetic {Esn = 3 x 10^^ erg) type lib SNe. As 
an alternative, J30l | suggested the importance of nuclei heavier than iron, which 
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are accelerated up to energies ~ 100 {Z\j = 92) larger than that of protons. 

We wish to stress however that the spectrum of CRs in this energy region is 
affected by the relative fluxes of Galactic and extragalactic CRs, and this makes 
the claim for excesses rather weak. The only point that appears rather strong 
in terms of the view of the origin of CRs illustrated in the papers listed above, 
as well as in the present paper, is that the spectrum of Galactic CRs is unlikely 
to extend to very high energies, thereby making the transition from Galactic 
to extragalactic CRs at the ankle rather poorly motivated. The two physical 
scenarios which are compatible with a transition in the energy region around 
lO'' — 10^ GeV are the mixed composition scenario [l|,0,y5lJl and the dip scenario 
111 . \l% y, Isl . The former model describes the extragalactic CR contribution as 



the superposition of different nuclei, so that the chemical composition in the 
transition region is mixed. The latter model is stunning for its simplicity: the 
extragalactic CR contribution is made of protons only and the transition from 
galactic to extragalactic CRs is completed at 10^^ eV. The transition occurs 
through a dip, produced by the onset of pair production, a feature which is very 
well defined in the CR spectrum. The chemical composition in the transition 
region changes suddenly [6[ from an iron dominated Galactic one to a proton 
dominated extragalactic one. 

In both scenarios the Galactic CR spectrum is cut off well below the ankle. 
The lowest energy part of the extragalactic CR spectrum is however affected, 
in both scenarios, by unknowns related to the propagation of extragalactic CRs 
in weak magnetic fields that might be present in the intergalactic medium, 
therefore the contribution of extragalactic sources to CRs in the energy region 
around a few x 10^ GeV remains poorly constrained. This uncertainty makes 
the need for an additional class of sources in the transition region somewhat 
weak. 



7. Discussion and conclusions 

The problem of explaining the origin of CRs as observed at the Earth is 
a complex one: it requires us to find the sources, describe the acceleration of 
nuclei (and in fact of electrons as well) in such sources, propagate all particles 
through the Galaxy taking into account diffusion and losses, and finally correct 
for local effects (such as the solar modulation at low energies). One can easily 
realize that this is a hard task, and yet one in which we have been successful in 
many respects, less in others. 

The most plausible sources of Galactic CRs remain SNRs, although a solid 
proof of the supernova paradigm satisfying all scientific standards has not been 
obtained as yet. The acceleration process is most likely diffusive shock ac- 
celeration in its non-linear version, which accounts for the dynamical reaction 
of accelerated particles on the shock and for the crucial phenomenon of self- 
generation of amplified magnetic fields in the acceleration region. In this re- 
spect we have been successful, in that several versions of the theory exist and 
the different formulations compare well with each other [see for instance ll9| . 
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Our semi-analytical approach has the advantage, when compared to all oth- 
ers, of simplicity and reduced computation time, which is the very reason why 
it is increasingly more adopted in complex hydrodynamical codes of the SNR 
evolution including the reaction of CRs. 

In the present paper we completed the theoretical framework previously put 
forward by our group with the inclusion of accelerated ions, which we find to 
contribute considerably to the shock modification. We show that the spectra 
of different chemicals observed at the Earth can be reproduced reasonably well 
and we also calculate the implications of acceleration of nuclei on the gamma- 
ray emission of an individual SNR. The all-particle spectrum of CRs is also 
reproduced in a satisfactory way. 

The one presented here is not the first attempt at describing acceleration of 
nuclear species in SNRs in the context of non-linear diffusive shock acceleration: 
[10[ published previous work on the topic using a numerical approach to the ac- 
celeration problem. Unfortunately a detailed comparison of our results with 
their work is made difficult by the fact that the procedure they used to include 
nuclei and the amplification of magnetic field are not fully illustrated: actually 
it is not even clear to us whether the nuclei were introduced in the non-linear 
chain or rather included as test particles in a shock structure mainly modified 
by CR protons. The spectra obtained in that paper for the individual nuclear 
species were very flat (and in fact concave) leading the authors to require a 
diffusion coefficient with a strong energy dependence, D{E) ex E^"^^. Such dif- 
fusion is however known to result in the breaking of the diffusive approximation 
well below the knee and consequently in excessive CR anisotropy, at odds with 
observations. 

More details were given in the work recently presented by Ptuskin et al. [38| . 
Their calculation consists again in a finite differences scheme for the solution 
of the coupled transport and fluid equations, as for |lO|. The authors discuss 
the effect of assuming a large velocity of the scattering centers, which leads to 
steepening of the resulting spectra (and to a reduction of the CR acceleration 
efficiency). As discussed in gl this very important issue, that was also discussed 
by Caprioli et al. |20| , is very poorly known: the velocity of the scattering centers 
can hardly be estimated reliably, and even more important the helicity of the 
waves is unknown. In fact, if the helicity is not chosen properly (but arbitrarily) 
the net effect may well be that of flattening the spectra even further. In the 
paper by Ptuskin et al. [38J, as well as in ours, it is assumed that an estimate 
of the speed of the waves is represented by the Alfven velocity evaluated in the 
amplified magnetic field at the shock. In J38l] the authors assume that the waves 
downstream of the shock move in a direction opposite to that of the shock. In the 
present paper we make the more conservative assumption that somehow waves 
are isotropized downstream, so that (vw) = 0. Both choices lead to spectra 
of accelerated particles which are steeper than the standard predictions of non- 
linear theory. In the absence of this very important ingredient it is not possible 
to explain CR spectra as observed at Earth. The main difference between our 
assumption and that made by Ptuskin et al. [38[ is that in the latter case 
the spectral slope is strongly dependent on the Alfvenic Mach number, which 
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leads the authors to suggest that some, yet unknown, self-regulation mechanism 
might_be at work to guarantee the appropriate level of field amplification. 

also made an attempt at taking into account the possibility that different 



classes of SNRs may contribute to the CR spectrum at Earth. However the 
relative abundances of chemicals injected in each class is the same and only 
the environmental conditions around the SNR (for instance gas density and 
temperature) are changed in the different cases. 

Although the introduction of these and other imaginable complications does 
not present any serious technical problem for the calculation procedure discussed 
here, we decided to focus on the main physical ingredients and avoid additional 
assumptions that, in our opinion, take away clarity from the results rather than 
making them more reliable. In our view, at the present stage of our knowledge, 
the assumptions that are required to carry out more sophisticated calculations 
are too many to lead to a scientifically trustable result. 

This leads us to a question that is of relevance for the investigation of the 
problem of CR origin from the theoretical point of view, namely "can we move 
beyond this point in a way that really adds to our knowledge of the problem?". 

There are clearly several open issues that are worth thinking about: (1) why 
are the observed spectra of protons and nuclei somewhat different? (2) what is 
the effect of the winds of the presupernova stars? (3) what is the best way to 
describe magnetic field amplification and/or complex topology of the magnetic 
field in the shock region? (4) how can we reliably describe the escape of CRs 
from a SNR? and (5) what is the Physics of injection of nuclei in a collisionless 
shock? 

It is possible (although far from guaranteed) that each one of these questions 
may be tackled from first principles, and if this ever happens, then the challenge 
will be putting together the conclusions in a convincing way to reproduce the CR 
spectra observed at Earth. Taking into account the wide variety of situations 
that are realized in Nature and the fact that in principle the CR spectra observed 
at Earth are the result of the convolution of so many different environmental 
situations, even this latter task would be formidable. At present, wc think it is 
difficult to envision going beyond the point already reached and still maintain 
the same standard of scientific credibility. On the other hand, the modeling of 
individual SNRs which leads to predict their appearance in different frequency 
ranges may provide several breakthroughs in the way SNRs accelerate nuclei and 
electrons, and this is probably the most promising avenue to follow for future 
progress. 

Acknowledgments 

We are grateful to the referee, Don Ellison, for his valuable comments. This 
work was partially supported by ASI through contract ASI-INAF 1/088/06/0. 



21 



References 

[1] Allard, D., Parizot, E., Olinto, A. V., Khan, E., Goriely, S., 2005, A&A, 
443, 29 



[6: 



[9 

[lo; 
[11 

[12: 
[13: 

[14 

[is: 

[16 

[17: 
[is: 

[19 

[2o: 

[21 
[22 
[23: 
[24: 



Allard, D., Parizot, E., Olinto, A. V., 2007, Astropart. Phys., 27, 61 

Allard, D., Olinto, A. V., Parizot, E., 2007, A&A, 473, 59 

Allard, D., Busca, N. G., Decerprit, G., Olinto, A. V., Parizot, E., 2008, 
JCAP, 10, 033 

Aloisio, R., Bcrezinsky, V., Blasi, P., Gazizov, A., Grigorieva, S., Hnatyk, 
B., 2007, APh, 27, 76 

Aloisio, R., Berczinsky, V., Blasi, P., Ostapchcnko, S., 2008, Phys. Rev. D, 
77, 5007 

Amato, E., Blasi, P., 2005, MNRAS Lett., 364, 76 

Amato, E., Blasi, P., 2006, MNRAS, 371, 1251 

BeU, A.R., 1978, MNRAS, 182, 147 

Berezhko, E. G., Volk, H.J., 2007, ApJ Letters., 661, 175 

Berezinsky, V., Gazizov, A., Grigorieva, S., 2005, Phys. Lett. B, 612, 147 

Berezinsky, V., Gazizov, A., Grigorieva, S., 2006, Phys. Rev. D., 74, 3005 

Blasi, P., Amato, E., Caprioh, D., 2007, MNRAS, 375, 1471 

Blasi, P., Gabici, S., Vannoni, G., 2005, MNRAS, 361, 907 

Bliimer, J., Engel, R., Horandel, J., 2009, PrPNP, 63, 293 

Caprioh, D., Blasi, P., Amato, E. & Vietri, M., 2008, ApJ Letters, 679, 139 

Caprioh, D., Blasi, P., Amato, E. & Vietri, M., 2009a, MNRAS, 395, 895 

Caprioh, D., Blasi, P., Amato, E., 2009b, MNRAS, 396, 2065 

Caprioh, D., Kang, H., Vladimirov, A., Jones, T.W., 2010, MNRAS, 417, 
1773 

Caprioh, D., Amato, E., Blasi, P., 2010a, APh., 33, 160 

Caprioh, D., Amato, E., Blasi, P., 2010b, APh., 33, 307 

Eichler, D., 1979, ApJ, 229, 419 

Ellison, D.C., Drury, L.O'C., Meyer, J-P., 1997, ApJ, 487, 197 

Ellison, D.C., Jones, F.C., Eichler, D., 1981, JGZG, 50, 110 



22 



[25] Ellison, D.C., Mobius, E., Paschmann, G., 1990, ApJ, 352, 376 

[26] Ellison, D.C., Patnaude, D.J., Slane, P., Raymond, J., ApJ, 712, 287 

[27] Higdon, J. C, Lingenfelter, R. E., 2005, ApJ., 628, 738 

[28] HiUas, A.M., 2005, J. Phys. G, 31, 95 

[29] Horandel, J. R., 2003, APh., 19, 193 

[30] Horandel, J. R., 2008, AdSpR, 41, 442 

[31] Horandel, J. R., Kalmykov, N. N., Tiniokhin, A. V., 2007, APh., 27, 119 

[32] Kaniae, T., Karlsson, N., Tsunefumi, M., Abe, T., Koi, T. 2006, ApJ, 647, 
692 

[33] Maestro, P., et al., 2010, Proceedings of the 31** ICRC, Lodz 2009 

[34] Malkov, M.A., Drury, L.O'C, 2001, Rep. Progr. Phys., 64, 429 

[35] Meyer, J-P., Drury, L.O'C, Ellison, D.C., 1997, ApJ, 487, 182 

[36] Meyer, J-R, Elhson, D.C., 1999, in LiBeB Cosmic Rays and Related X- 
and Gamma-rays, ASP Conference Series, Vol. 171, eds., Ramaty, Vangioni- 
Flam, Cass & Olive, p. 187 

[37] Ptuskin, V.S., Zirakashvih, V.N., 2005, A&A, 429, 755 

[38] Ptuskin, V.S., Zirakashvih, V.N., Seo, E.S., 2010, ApJ, 718, 31 

[39] Skilling, J., 1975, MNRAS, 172, 557 

[40] Truelove, J. K. and Mc Kee, C. F., 1999, Apj Supplement Series, 120, 299 

[41] Wiedenbeck, M.E. et al., 2007, Space Sci. Rev., 130, 415 



23 



